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Gravitino Dark Matter represents a compelling scenario in Supersymmetry, which brings together 
a variety of data from cosmology and collider physics. We discuss the constraints obtained from the 
LHC on supersymmetric models with gravitino dark matter and neutralino NLSP, which is the case 
most difficult to disentangle at colliders from a neutralino LSP forming DM. The phenomenological 
SUSY model with 19+1 free parameters is adopted. Results are obtained from broad scans of the 
phase space of these uncorrelated parameters. The relation between gravitino mass, gluino mass 
and reheating temperature as well as the derived constraints on these parameters are discussed 
in detail. This relation offers a unique opportunity to place stringent bounds on the cosmological 
model, within the gravitino dark matter scenario, from the results of the LHC searches in Run-2 
and the planned High-Luminosity upgrade. 

PACS numbers: 11.30.Pb, 14.80.Ly, 14.80.Nb, 95.35.+d 


I. INTRODUCTION 

With the discovery of the Higgs boson at the LHC [L, 
[2], collider particle physics moves on to a new phase, 
where searches for new phenomena will be at the center 
of the physics program. Dark matter (DM) plays a spe¬ 
cial role in connecting these searches to cosmology. Not 
only DM appears to be the most solid indication of new 
physics obtained so far, but the precision of the deter- 
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mination of its relic density in the Universe, through the 
study of the cosmic microwave background (CMB) ||3|, 
and the increasing sensitivity of the searches for DM scat¬ 
tering in underground experiments ms] represent power¬ 
ful constraints for most models of new physics that offer 
a DM candidate. 

Supersymmetry (SUSY) is one of the most promising 
and best motivated extensions of the Standard Model of 
particle physics (SM), providing at the same time a so¬ 
lution to the hierarchy problem within a weakly coupled 
theory, radiative electroweak symmetry breaking with a 
light Higgs scalar with properties that, in the decoupling 
limit, come to resemble those of the Standard Model, and 
also naturally includes several viable dark matter candi¬ 
dates within its extensive particle spectrum. 

Searches conducted at the LHC on the 7+8 TeV Run-1 



2 


data have already set some significant constraints on low 
energy SUSY, in particular on the mass of the coloured 
superpartners and on the general mass scale in the case of 
highly constrained models, such as the CMSSM and the 
NUHM models [7-9]. However, many scenarios with su¬ 
persymmetric particles within the reach of the LHC runs 
at 13-14 TeV are totally unprobed at present, in partic¬ 
ular when we consider more general SUSY models EH - 
[13] such as the “phenomenological MSSM” (pMSSM), 
without any assumption on the SUSY breaking media¬ 
tion mechanism and on unification at high scale. This 
model gives the most general framework in the minimal 
SUSY extension of the SM (MSSM) for studying super- 
symmetric signals and constraints from colliders. It is 
being adopted by an increasing number of phenomeno¬ 
logical mm and experimental OEU studies. 

The particle most studied in conjunction with SUSY 
dark matter is the lightest neutralino, y 9 , which pro¬ 
vides us with a viable realisation of the WIMP mecha¬ 
nism and implies DM signals not only at colliders, but 
also in direct and indirect detection. However, another 
very well-motivated SUSY candidate for DM is the grav- 
itino, G, the superpartner of the graviton which couples 
with Planck-suppressed strength to the rest of the model, 
making it very difficult to detect in direct or indirect DM 
experiments, in SUSY models with conserved R-parity. 
Even a small R-parity breaking opens up the possibil¬ 
ity of gravitino decaying DM candidate relaxing part of 
the cosmological constraints [28] . Colliders offer a unique 
way to probe the DM sector for scenarios with gravitino 
dark matter. Therefore, dedicated searches should be 
actively pursued to fully exploit the LHC reach in these 
scenarios. 

In general, assuming that the gravitino is the Lightest 
SUSY Particle (LSP) and forms DM greatly relaxes the 
DM constraints on the SUSY spectrum, since its abun¬ 
dance is not produced through the WIMP mechanism 
and does not only depend on the SUSY spectrum, but 
also on the reheating temperature in the Universe after 
inflation naisni- Consequently, the G remains a good 
candidate for DM even if the spectrum of the superpart¬ 
ners turns out to be much heavier than expected. Nev¬ 
ertheless, we will show that the correlation between the 
gluino (and in general gaugino) mass and gravitino pro¬ 
duction ensures that the LHC will be able to probe the 
gravitino production mechanism for high reheating tem¬ 
peratures above 10 9 GeV and that the High Luminosity 
LHC program (HL-LHC) may play a crucial role for these 
tests. 

In this paper we present a study of the present and fu¬ 
ture constraints on SUSY models with gravitino LSP re¬ 
sponsible for DM and neutralino Next-to-Lightest SUSY 
Particle (NLSP). These models are not yet strongly 
constrained since the current LHC limits on weakly- 
produced SUSY particles are limited in coverage and our 
scenarios of interest not easy to disentangle from those 
of neutralino LSP. Indeed, here we consider the parame¬ 
ter space where the neutralino is stable on collider time 


scales and gives rise to the usual missing transverse en¬ 
ergy (MET) signatures at the LHC. After considering the 
present bounds on the model and how strongly the cos¬ 
mological and astrophysical constraints limit the param¬ 
eter space, we assess the capability of the forthcoming 
LHC runs to test the gravitino production mechanism 
and to tell the gravitino LSP from the neutralino LSP 
solutions. 

Compared to previous studies, this analysis imple¬ 
ments the latest collider constraints, including those from 
mono jet searches which complement the other searches 
in the regions of the parameter space with degenerate 
SUSY masses, and we discuss the combination of LHC 
and other data in the context of the gravitino models. We 
adopt the pMSSM as a generic MSSM model that, con¬ 
trary to more constrained models used in the past [3Th37] 
for similar studies, does not imply relations between the 
masses of the different SUSY particles. This opens up 
new region of the parameter space with interesting phe¬ 
nomenology and search opportunities at the LHC which 
were not available in the constrained MSSM (CMSSM), 
as we demonstrate in this study. The concurrent analy¬ 
sis of the LHC data and the DM direct detection experi¬ 
mental results will be key to identify a neutralino NLSP 
signal without observation of the decay or to exclude a 
large fraction of the model parameter space. The combi¬ 
nation of the LHC sensitivity on the gluino mass and the 
requirement of a large enough reheating temperature af¬ 
ter inflation for thermal leptogenesis [38H4T] is specific to 
the study performed in the pMSSM, restricts the viable 
MSSM parameter space and highlights the capability of 
the LHC to test these cosmological scenarios. 

This paper is organised as follows: after discussing the 
generalities of the gravitino LSP and DM scenario to¬ 
gether with the cosmological bounds from nucleosynthe¬ 
sis in chapter II, we describe our scan strategy for the 
pMSSM parameter space and the collider and low-energy 
constraints in chapter III. Chapter IV presents the results 
of our analysis, while in chapter V we give our conclu¬ 
sions. 


II. GRAVITINO DARK MATTER 

The gravitino is the superpartner of the graviton in 
models with local supersymmetry 02- It is the particle 
most directly related to the effect of SUSY breaking and 
obtains its mass via the SuperHiggs mechanism [43] as 
the gauge-fermion of SUSY from any existing source of 
breaking. The hierarchy between the masses of the super¬ 
partners is determined by the particular mediation mech¬ 
anism, and it can result in a gravitino which is heavier, 
lighter or even much lighter than the other superpartners. 
If the gravitino is the lightest SUSY state and therefore 
absolutely (or sufficiently) stable, it represents a viable 
candidate for explaining dark matter. The gravitino cou¬ 
plings are set by supergravity and by the MSSM parame¬ 
ters and are inversely proportional to the reduced Planck 
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mass, taken here at its standard value of Mp = 2.4 x 10 18 
GeV. The gravitino mass is the only additional parame¬ 
ter needed to describe the gravitino and its interactions. 
Starting with the phenomenological MSSM with 19 pa¬ 
rameters set at the electroweak scale, we will add Mq as 
independent parameter, bringing the total to 20 . 

In this paper, we shall focus our analysis on the part of 
the parameter space where the gravitino is the LSP and 
a cold DM candidate. We do not consider a very light 
gravitino, with mass in the eV to few keV range, since it 
is at odds with present cosmological data which strongly 
disfavour hot or warm Dark Matter and may only con¬ 
stitute a subdominant DM component [44] [47 . In this 
region of parameters, another (cold) DM candidate is 
needed and no other supersymmetric particle can fulfill 
this role. Moreover, a very light gravitino can be a warm 
thermal relic only in presence of entropy production with 
a sufficient dilution factor of the gravitino freeze-out den¬ 
sity [48'. Still, this scenario has been widely studied, 
in particularly in relation to the gauge-mediated SUSY 
breaking (GMSB) mechanism 49, 50] . It gives very dis¬ 
tinctive signatures at colliders since the NLSP decays 
promptly leading to photons/jet + MET, which are ac¬ 
tively searched by the ATLAS and CMS experiments at 
the LHC [5ffl54] . In the scenarios we explore here, the 
Next-to-Lightest SUSY Particle (NLSP) often mimicks a 
stable LSP, even if the NLSP can be in principle charged 
under the SM group and therefore cannot be stable on 
cosmological scales. We assume the NLSP to be the light¬ 
est neutralino, which is the case most difficult to dis¬ 
entangle at colliders from the scenario of neutralino LSP 
and Dark Matter particle. 

Although very weakly coupled, gravitinos can be pro¬ 
duced in substantial numbers in the primordial thermal 
bath from scattering of SUSY particles, in particular the 
colored or electroweakly charged states. Usually the in¬ 
teractions involving gluinos or gauginos gives the domi¬ 
nant contribution since they are some of the most abun¬ 
dant and strongly coupled particles and interact with 
gravitinos by an effective dimension-5 operator, leading 
to a yield proportional to the highest equilibrium temper¬ 
ature of the thermal bath [29| 130 ] . For a light gravitino, 
this contribution is dominated by the Goldstino compo¬ 
nent of the gravitino and given by 


nV? h 2 = 0.83 


TrH / w 3/2 \ 1 s ( Mj 

10 9 GeV V 1 GeV/ RH ’ \300GeV 


( 1 ) 

where 7 iiTpp) are numerical factors including the depen¬ 
dence on the gauge coupling for each SM gauge group and 
the evolution of the gauge couplings and gaugino masses 
Mi to the scale of the reheating temperature Tpp. To 
do so, we used here only one-loop RGEs to capture the 
main effect and we define, following m, 


li(T) = yi 


gf (T) 


log 


gf(Mi) [gi(T)\ 


ki 


( 2 ) 


with (yi,ki) equal to (4.276,1.271), (1.604,1.312), 


(0.653,1.266) for the SM gauge groups SU( 3), SU( 2 ) and 
1/(1), respectively. It is clear from these factors that the 
gluinos often give the largest contribution, due to their 
multiplicity and their larger gauge coupling. The ex¬ 
pression given in eq. 0 is valid only in the perturbative 
regime, and becomes less reliable at low reheat temper¬ 
atures, where g 3 becomes large. Note though, that for 
low Trh , in the region of supersymmetric masses we will 
explore, gaugino scatterings do not provide the dominant 
contribution to the gravitino relic density and therefore 
an order one determination of this part is sufficient for 
our purposes. 

Processes involving scalar SUSY particles, like squarks 
and sleptons, include also dimension-4 type couplings, 
which generate a yield independent of the thermal bath 
temperature, but those contributions are often negligi¬ 
ble, unless a strong hierarchy between scalar and gaug¬ 
ino masses is assumed [55 , 5b. The contribution from the 
NLSP decay after freeze-out can be quite substantial, es¬ 
pecially for the case of a Bino-like neutralino NLSP. Such 
contribution is given by [57H59] 


n^h^(^y x h\ ( 3 ) 

where M x is the mass of the lightest neutralino and h 2 
its present energy density, if the neutralino would be sta¬ 
ble, set by the freeze-out process. 

In determining the relic gravitino density we consider 
the thermal contribution from gaugino scattering given 
in eq. 0 as well as that arising from neutralino decay out 
of equilibrium in eq. ©• For any value of the pMSSM 
parameters and gravitino mass, a specific value of the re¬ 
heat temperature allows us to match the observed DM 
abundance, as long as the contribution from NLSP decay 
is smaller than the observed DM density. In the pMSSM 
parameter region where the neutralino relic density ex¬ 
ceeds the DM energy density, an upper bound on the 
gravitino mass also appears. This bound may be stronger 
than our request that the gravitino is the LSP. In general, 
there is a maximal value of the reheating temperature al¬ 
lowed at any parameter point, and possibly also a max¬ 
imal gravitino mass, in order to avoid gravitino density 
overclosure. 

The total gravitino density can be written in a simpli¬ 
fied form as 


Vt^h 2 = A 


Trh 

m 3 / 2 


+ B m 3 /2 


( 4 ) 


where the parameters A and B are functions of the 
pMSSM parameters, as can be read from eqs. 0 ,@. A 
depends mostly on the gaugino masses Mi, M 2 and M 3 , 
while B contains the neutralino relic density and depends 
directly on the neutralino composition. Requiring these 
two mechanisms to give rise to the full dark matter en¬ 
ergy density can be recast in a simple quadratic equation 
for m 3/2 : 

B ml/2 ~ ^ CDMh 2 TO3/2 + A Trh = 0 ( 5 ) 
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which has two real roots for ra 3 / 2 only when 


A. Cosmological constraints 


Trh < 


(flcDM ^ 2 ) 2 
4 AB 


(6) 


and this value provides the maximal possible reheat¬ 
ing temperature for any gravitino mass. Note that the 
r.h.s. of eq. § depends only on the pMSSM param¬ 
eters through the A and B coefficients. The maximal 
reheating temperature is realized when the CDM density 
is provided equally by the two production mechanisms 
and the gravitino mass is fixed to the value 


%/2 = 


^CDAlh 2 
2 B 


^CDM h 2 M x 

Q x h 2 Y 


(7) 


If Q x h 2 < ftcDMh 2 / 2 , this particular value of ra 3 / 2 is ex¬ 
cluded by the requirement ra 3 / 2 < M x . For values of the 
gravitino mass lower than in eq. 0/ the bound on Trh 
gets stronger. Other cosmological constraints restrict the 
region of large ra 3 / 2 and in general in any point of the 
pMSSM parameter space set a stronger bound on Trh 
compared to that given by eq. (|6|. This constraint is 
given by 


Trh < {p'CDMh 2 — Bm 3 / 2 ) , (8) 

and it depends on both the gravitino mass and the 
pMSSM parameters. 

The reheating temperature of the Universe is an im¬ 
portant parameter in the cosmological evolution and de¬ 
pends on both the particular inflationary model and the 
reheating process after inflation. It strongly influences 
the possible mechanisms of baryogenesis, needed to gen¬ 
erate the present baryon asymmetry starting from a sym¬ 
metric initial state. In particular, the thermal leptogen- 
esis mechanism m relies on the presence of heavy right- 
handed Majorana neutrinos in the thermal bath. These 
decay out-of-equilibrium giving rise to a non-vanishing 
lepton and then a baryon asymmetry. 

Thermal leptogenesis can produce the observed baryon 
number only if the reheating temperature is sufficiently 
high, above approximately 2 x 10 9 GeV [38l[40]. With 
a mild tuning in the seesaw formula and exploiting fla¬ 
vor effects, this bound can be relaxed by an order of 
magnitude m eb and even more so in case of reso¬ 
nant CP violation, due to nearly degenerate RH neutrino 
masses [62, 63] . While these bounds are not intrinsic to 
the gravitino DM hypothesis, they are a desirable addi¬ 
tion to Big Bang cosmology. We will discuss how the 
requirement of high Trh can be used in these specific 
scenario to place strong constraints on the SUSY spec¬ 
trum, when it is compared to the reheating temperature 
needed to generate the correct gravitino DM abundance. 
Indeed, as can be inferred from eq. 0, light degener¬ 
ate gauginos reduce gravitino production and allow for 
the highest possible Trh around 10 9 GeV [64, 65 . This 
makes possible to interpret collider data as tests of the 
simplest thermal leptogenesis hypothesis. 


The requirement of a consistent cosmology sets impor¬ 
tant constraints on models with a very weakly-interacting 
particle like the gravitino. Indeed, even if the gravitino 
is not the LSP and therefore unstable, it may decay very 
late in the cosmological history leading to the disrup¬ 
tion of Big Bang Nucleosynthesis (BBN), the so-called 
”gravitino problem” [66H68] . It is interesting to point 
out that these scenarios may be seen as an asset, since 
these late decays might solve another outstanding prob¬ 
lem, that of the amount of cosmic Lithium formed af¬ 
ter the Big Bang [69]. In our scenario the gravitino 
is the LSP and so the long-lived particle becomes the 
neutralino NLSP. Due to the Planck scale suppression, 
the NLSP lifetime can easily exceed 1 s and the de¬ 
cay of the out-of-equilibrium NLSPs again takes place 
during BBN, thus affecting the abundance of light ele¬ 
ments EH- For a neutral long-lived particle like the neu- 
tralino, the main processes affecting BBN are photo- or 
hadro-dissociation, which depend on the number density 
of the decaying particle and its branching fraction into 
photon or hadrons. For short lifetimes the constraints 
from hadro-dissociation are stronger than those given by 
photo-dissociation, which takes over at r > 10 4 — 10 6 
s [ZH [72]. The constraints for neutral relic with differ¬ 
ent values of the hadronic branching ratio and decaying 
particle mass were computed by running a BBN code 
including the neutral particle decay channel into quark 
and antiquark pairs in [72] . Here, we use those generic 
constraints and translate them to the special case of the 
neutralino, similarly as what was done in pang. In or¬ 
der to do that we need to compute the decaying branch¬ 
ing fractions and the freeze-out density of the NLSP for 
every point of the pMSSM parameter space. 

The decay rates of a general neutralino NLSP into a 
gravitino LSP were calculated analytically including all 
possible three-body decays in J73j. The neutralino decays 
dominantly into a gravitino and a gauge boson. The 
hadronic channel into gravitino with a quark and anti¬ 
quark pair arises most of the time from an intermediate 
gauge or Higgs boson, as long as the scalar quarks are not 
too light. So the most important hadronic channel for the 
neutralino consists of a three-body decay, where part of 
the energy is carried away by the inert gravitino, instead 
of the simple two-body decay into quark and antiquark, 
which was assumed in m for the neutral relic. There¬ 
fore the quark-pair momentum distribution in neutralino 
decay is softer than in a two-body decay and we there¬ 
fore regard the limits given in m as possibly too strong 
for our scenario. Indeed, if the intermediate gauge boson 
is on-shell before decaying into a quark-antiquark pair, 
the number of hadronic particles produced is expected 
to be independent of the gauge boson energy or momen¬ 
tum. Moreover, if the hadrons produced thermalise be¬ 
fore interacting with light elements, any dependence on 
the initial spectrum or the mass of the mother particle 
is washed out. In this spirit, the work of ref. m uses 
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the bounds of Jedamzik for a decaying particle mass of 
100 GeV, disregarding the dependence on this mass and 
correcting the electromagnetic bounds by subtracting the 
energy carried away by the gravitino. In this study we 
implement both sets of constraints on our scans: first 
the limits given by m as a function of the hadronic 
branching ratio and decaying particle mass, and then 
the slightly weaker limits of ref. m- Indeed, since in 
this analysis we include also regions of the parameter 
space where the Z or Higgs boson cannot be on shell, 
i.e. the region where Mi — ra 3 / 2 < M^, M^, we consider 
the latter bounds as too conservative in some part of the 
parameter space. The comparison between the two sets 
of bounds, with one appearing to be too strong and the 
other too weak in part of the parameter space, provides 
us with an estimation on the sensitivity of the results on 
the precise implementation of these constraints. 

In general, BBN bounds exclude regions where the den¬ 
sity of the decaying particle as a function of its lifetime is 
too large. Indeed in the early stages of BBN, correspond¬ 
ing to lifetimes below 10 2 s, the constraints are very weak, 
or even non-existent, while after the production of light 
elements, like Deuterium, Helium or Lithium, the effect 
of hadro-dissociation becomes stronger. The lifetime of 
the NLSP in our scenario is always related to the grav¬ 
itino mass, since for a light gravitino its couplings are 
proportional to l/m 3 / 2 . These bounds can be translated 
in general into an upper bound on the gravitino mass, as 
a function of the NLSP density at freeze-out. For a bino- 
like neutralino NLSP lifetime, above the Z threshold, we 
have 



1 


- X?-n/^b3/2 




487 xMp rn*/2 


= (57s) 


-l 


f M X V, 

( m 3/2 W 

V 1 TeV J 

U0 GeV fin 


(9) 

For gravitino masses larger than about 1 GeV, the NLSP 
typically decays during BBN and affects its predictions. 
It is possible to shorten the NLSP lifetime by varying the 
value of the gravitino mass, independently of the decay¬ 
ing neutralino mass or composition. 

We use the results of m to implement the hadronic 
constraints on the model, deriving the branching frac¬ 
tion for the neutralino decays into gravitino and quark- 
antiquark pair for each point of the pMSSM parameter 
space scan. This branching fraction is shown as a func¬ 
tion of the neutralino mass for our pMSSM scan points 
in Figure [l] The nature of the neutralino is highlighted, 
with the bino-, wino- and higgsino-like components de¬ 
fined as the points with the corresponding entry squared 
in the neutralino mixing matrix exceeding 0.75. The 
minimal hadronic branching fraction for fixed neutralino 
mass, is obtained for the neutralino composition of a pure 
photino, since in this case the contributions through in¬ 
termediate Z and Higgses are absent. 

We obtain the neutralino branching fraction into 
hadrons as follows. First, we estimate the neutralino 
lifetime considering the decay into the gravitino and a 
photon, which is always open, keeping all kinematical 


FIG. 1: Neutralino hadronic decay branching fraction as a 
function of the neutralino mass for the accepted pMSSM scan 
points, with the nature of the neutralino having the largest 
contribution in each region highlighted. 


factors in this channel. Then, we compute all the other 
decay widths, e.g. into Z and h, by setting the gravitino 
mass to zero in the kinematical factors. Since the photon 
or another open decay channel dominates, when thresh¬ 
olds are encountered, this approximation holds below the 
per-mille level in all points and gives a rather conserva¬ 
tive estimate of the hadronic branching fraction. Next, 
e compute the decay rate into the hadronic channel tak- 
g into account the neutralino composition and the pres¬ 
ence of off-shell or on-shell intermediate bosons. If the 
intermediate particles are on-shell, the result is well ap¬ 
proximated by the decay rate into that particle times the 
corresponding branching fraction. But our results extend 
also into the regions with off-shell intermediate states. In 
this case, we integrate the branching fraction for energies 
of the quark-antiquark pairs above their mass threshold. 
As discussed in m this gives, below the Z threshold, 
a result three times larger than taking a threshold of 
2 GeV, corresponding to the production in the decay of 
hadrons instead of mesons. At leading order the branch¬ 
ing fraction is independent of the gravitino mass, since all 
rates are proportional to the same factor l/m 2 ^ 2 . Only 
for large gravitino masses near thresholds does this pa¬ 
rameter play an important role also in the phase-space 
factors. 

Photo-dissociation constraints are very similar in most 
of the parameter space and mostly independent on the 
specific values of the purely electromagnetic branching 
ratios, since the decay into the gravitino and a pho¬ 
ton has always a branching fraction larger than 30% 
and the hadronic bounds in ref. [72] include the photo¬ 
dissociation effects from charged hadrons. 

The calculation of the cosmological bounds require to 
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implement the number density of the neutralino NLSP 
before decay. This is evaluated by considering the freeze- 
out process of the neutralino and solving the corre¬ 
sponding Boltzmann equation using the Super I so Relic 
code m- The neutralino NLSP density is determined 
by the 19 pMSSM parameters and it is independent of 
the gravitino mass. Therefore we can find iteratively the 
largest NLSP lifetime and thus the largest gravitino mass 
compatible with BBN and the NLSP number density for 
a set of pMSSM parameters, as discussed in the next sec¬ 
tion. The maximal allowed value for the gravitino mass 
at any chosen point in the pMSSM parameter space is 
usually smaller than the general upper bound given in 
eq. 0 and therefore we can apply eq. ^ to obtain a 
tighter bound on the reheating temperature Trh- The 
distribution of the neutralino relic density as a function 
of its lifetime for the pMSSM scan points is given is Fig¬ 
ure [2] The dependence of the relic density of the neu- 
tralino from its composition clearly appears: a bino neu- 
tralino has naturally a larger abundance and therefore as 
an NLSP must have a shorter lifetime than a Wino and 
Higgsino neutralino in order to not disrupt BBN. We ex¬ 
pect that wino/higgsino-like NLSPs can allow for a larger 
maximal gravitino mass, and so also a larger reheating 
temperature than a bino NLSP. This can be verified on 
our scans by studying the distribution of the reheating 
temperature as a function of the NLSP mass and select¬ 
ing the NLSP nature. Figure [3] shows the results ob¬ 
tained constrasting the case of the pMSSM scan to those 
of the cMSSM. In the pMSSM parameter space, points 
with wino- and higgsino-like NLSP accomodate values of 
T rh in the range 10 6 - 5xl0 9 GeV. On the contrary, the 
bino-like population of the cMSSM points does not offer 
solutions with T RR values in excess of 10 7 GeV. 

For neutralino lifetimes longer than 10 10 s, constraints 
from CMB distortion also play a role and start to be 
more stringent than the BBN constraints, see e.g. m- 
Since we can obtain such extremely large lifetime values 
only for very few fine-tuned points of our parameter scan, 
as in the degenerate gravitino-NLSP scenario m , we 
implement an upper cut at r = 10 11 s. 


III. ANALYSIS 
A. MSSM Scans 

The phenomenological MSSM (pMSSM) is the most 
general MSSM scenario with R-parity, CP conservation 
and minimal flavour violation m and its properties are 
defined by 19 independent parameters. This analysis is 
based on a flat scan where these parameters are varied 
within the ranges given in Table [TJ The gravitino mass 
is added as an additional parameter. The ranges of the 
SUSY particle masses have been made broader compared 
to previous studies [16] [18] [19] to achieve a better assess¬ 
ment of the capability of the LHC at 13-14 TeV, which 
will probe the masses of strongly interacting SUSY par- 




FIG. 2: Neutralino relic density as a function of the neutralino 
lifetime for the accepted pMSSM scan points, with the nature 
of the neutralino having the largest contribution in each re¬ 
gion highlighted. The BBN constraints obtained following 
ref. [72] (top) and ref. 19] (bottom) are shown by the con¬ 
tinuous lines corresponding to limits for neutralino hadronic 
branching fractions of 100%, 10%, 1% and 0% from bottom 
to top. 

tides up to 2-3 TeV. For this analysis a total of 7xl0 7 
pMSSM points have been generated and analysed. 

In our analysis, the SUSY mass spectra and couplings 
are generated from the input pMSSM parameters with 
SOFTSUSY 3.3.3 Only points having the lightest 

neutralino, y?, to be the next lightest SUSY particle to 
the gravitino and the lightest Higgs mass in the range 
122 - 128 GeV are accepted. Particle decay widths are 
calculated using HDECAY 5 m and SDECAY m • 

In order to contrast the gravitino phenomenology in 
the parameter space of the pMSSM to that of constrained 
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FIG. 3: Reheating temperature as a function of the NLSP 
neutralino mass for the accepted pMSSM (top) and cMSSM 
(bottom) scan points, with the nature of the neutralino having 
the largest contribution in each region highlighted. 


MSSM models, we also perform scans of the constrained 
MSSM (CMSSM) gal |8l] where we vary the Mo, M]y 2 , 
tan f3 and Aq parameters in the ranges indicated in Ta¬ 
ble HI 

For each accepted point, we compute the flavour ob¬ 
servables and the muon anomalous magnetic moment 
with Superlso v3.3 [82l [83], and the neutralino relic 
density with Superlso Relic [74] . where the gravitino 
related calculations have also been implemented. To 
allow for comparisons between the gravitino LSP and 
the neutralino LSP scenarios, we compute the neu¬ 
tralino direct detection scattering cross sections using the 
MicrOMEGAs code [84] . 

The constraints from the LHC SUSY searches at 7+8 


Parameter 

Range 

tan [3 

[1, 60] 

Ma 

[50, 5000] 

Mi 

[-5000, 5000] 

m 2 

[-5000, 5000] 

M 3 

[0, 5000] 

A d = A s = A b 

[-15000, 15000] 

II 

O 

II 

[-15000, 15000] 

A e = A„ = A T 

[-15000, 15000] 


[-5000, 5000] 

Me L — M^ l 

[0, 5000] 

Me R — Mfx R 

[0, 5000] 

Mf L 

[0, 5000] 

Mf R 

[0, 5000] 

Mq 1L Mq 2L 

[0, 5000] 

Mq 3L 

[0, 5000] 

Mu R — m~ cr 

[0, 5000] 

Mr 

tR 

[0, 5000] 

II 

sT 

[0, 5000] 

Mr 

_ks_ 

[0, 5000] 


TABLE I: Range of the pMSSM parameters adopted in the 
scans (in GeV when applicable). 


Parameter 

Range 

tan f3 

[2, 68] 

Mo 

[0, 10000] 

Mi/2 

[0, 10000] 

Aq 

[-10000, 10000] 


TABLE II: Range of the cMSSM parameters adopted in the 
scans (in GeV when applicable). Both signs of /i have been 
considered. 


and their projection to 14 TeV are obtained through the 
analysis of inclusive SUSY events generated by PYTHIA 
8.150 [85] with the CTEQ6L1 parton distribution func¬ 
tions [SB. SUSY particle production accompanied by a 
hard jet, relevant to mo no jet searches, is simulated using 
MadGraph sm followed by PYTHIA for hadronisation, as 
described in [24]. The physics objects of the signal events 
are obtained with a parametric simulation of the LHC 
detector response performed with Delphes 3.0 [88]. 


For each pMSSM point, we infer the maximal grav¬ 
itino mass consistent with the BBN constraints and lead¬ 
ing to a relic density QqJi 2 ^ 0.11, in agreement with 
the Planck CMB data [3+ First, we derive the maximal 
neutralino lifetime using the BBN limits given in sec¬ 
tion II A| from the computed neutralino relic density and 
hadronic branching ratio. The corresponding gravitino 
mass is then obtained through an iterative procedure 
from the maximal NLSP lifetime. Once the maximal 
gravitino mass is obtained, we randomly scan over the 
gravitino mass values within a range comprised between 
1/1000 and 10 times the maximal gravitino mass and 
keep ten of such gravitino mass values for each pMSSM 
point. If the maximal gravitino mass cannot be obtained, 














































Flavour 

2.63 x 10“ 4 < BR(B -> X sl ) < 4.23 x 10“ 4 [M] 

1.3 x 10“ 9 < BR(B S -> p + /Ountag < 4 - 5 x 10 “ 9 EDHSU 
0.33 x 10“ 4 < BR (B u -> tv) < 1.95 x 10“ 4 [H 02] 

4.7 x 10“ 2 < BR (D s -s- tv) < 6.1 x 10“ 2 [H§][M] 
0.985 < 23 < 1-013 [96] 

G 9 - 2)„ 

-2.4 x 10" 9 < Sa^ < 5.0 x 10“ y [1302] 


TABLE III: Summary of the constraints from flavour physics 
and muon anomalous magnetic moment. 


we vary randomly the gravitino mass logarithmically be¬ 
tween 0.1 MeV and the neutralino mass. We then com¬ 
pute the gravitino relic density from NLSP decay as given 
in eq. and check if it is smaller than or equal to the 
DM relic density. Since the neutralino relic density is 
diluted by the ratio of the gravitino to the neutralino 
mass and the BBN bounds are stronger for larger NLSP 
density, this is always the case. To retrieve a relic den¬ 
sity compatible with the cosmological observations, we 
compute by inverting eq. ([l]) the reheating temperature, 
Trh , required to increase the gravitino abundance up to 
the value corresponding to the Planck observations. 

In addition, we include a set of points for which the 
lightest neutralino is a photino, i.e. a mixture of bino and 
wino with the dominant decay mode to a gravitino and 
a photon. These states can be obtained at tree level for 
Mi = M 2 and |Mi| <C |/i|. However, the neutralino mass 
matrix receives higher order corrections, which modify 
this condition. Therefore, we use a specific scan where 
we first impose Mi = M 2 , and then iteratively adjust the 
value of M 2 to retrieve the pure photino solution. 


B. Collider and low-energy constraints 

The analysis of the viable MSSM parameter space re¬ 
quires the implementation of the available constraints 
from different sectors. 

First the constraints from low energy observables are 
imposed to the generated points. Table m summarises 
the most recent values for the low energy 95% C.L. lim¬ 
its, obtained by including experimental and theoretical 
uncertainties. In particular, the branching fraction of 
the decay B s -+ /i~ is sensitive to scalar and pseu¬ 

doscalar operators, and provides us with a bound on the 
tan/3 and Ma parameters [98] . The value of the inclu¬ 
sive BR (B -+ X s +) is also sensitive to charged Higgs 
and charginos/stops contributions at loop level, thus con¬ 
straining tan/3, Ma, as well as the stop and chargino 
masses. These constraints remove ~ 12% of the accepted 
pMSSM points in our scans. 

With ~25 fb -1 of statistics collected by both the AT¬ 
LAS and CMS experiments at 7+8 TeV during the LHC 
Run-1, a vast array of searches for the production and 
decay of new particles has been performed. In particu¬ 
lar, the analysis of channels with missing transverse en¬ 
ergy (MET), including jets + MET, leptons + MET and 


Channel 

Experiment 

Sensitivity 

Ref. 

jets + MET 

ATLAS 

9, 9 

TST 

2 1 + jets + MET 

ATLAS 

t,q 

.100. 

1 £ + b-jet(s) + MET 

ATLAS 

t 

101 

6-jets + MET 

ATLAS 

7, b 

102 

2 £ + MET 

ATLAS 

x u x ± 

.103 

3 £ + MET 

ATLAS 

x° x ± 

.104 

1 £ + bb + MET 

ATLAS 

x° x ± 

105 

mono jet + MET 

ATLAS 

XX, qq 

.106 

mono jet + MET 

CMS 

xx, qq 

107 

H/A -+ TT 

CMS 

H, A 

108| 


TABLE IV: Summary of the analyses used to assess the ob¬ 
servability of the pMSSM points by the LHC SUSY searches. 


monojets, have excluded a sizeable fraction of the MSSM 
parameter space. With the start of Run-2 at 13 TeV and 
the perspectives of increased luminosity, the LHC will put 
SUSY scenarios through a crucial, although possibly not 
yet definitive, test. In this study, we consider the bounds 
obtained from the analyses summarised in Table |IV[ Al¬ 
though these are not exhaustive of the signal topologies 
investigated by the ATLAS and CMS searches, they cover 
the SUSY signatures with the highest sensitivity and are 
largely uncorrelated. For each accepted pMSSM point, 
we simulate an inclusive SUSY event sample and perform 
a parametric simulation for the event reconstruction, as 
discussed above. Signal selection cuts corresponding to 
each of the analyses are applied to these simulated signal 
events. The number of SM background events in the sig¬ 
nal regions are taken from the estimates reported by the 
experiments. The 95% confidence level (C.L.) exclusion 
in presence of background only is determined using the 
CLs method 

These results are projected to 14 TeV for 300 and 
3000 fb -1 of integrated luminosity, by generating events 
at 14 TeV and rescaling the 8 TeV backgrounds by the 
corresponding increase in cross section and signal cut ac¬ 
ceptance at the higher energy. The constraints derived 
from these searches exclude a significant fraction of the 
gravitino pMSSM points, from 22% at 8 TeV to 75% 
and 85% at 14 TeV with 300 and 3000 fb -1 , respectively. 
The gluino mass is of particular importance in relation 
to the gravitino relic density and the reheating tempera¬ 
ture. Figure [4] shows the distribution of the gluino mass 
for the points not excluded by the searches for the differ¬ 
ent energies and data sets considered here. The gluino 
mass values at which more than 95% of the scan points 
below that mass are rejected by the LHC searches, or will 
be in case of a negative result, is Mg = 840, 1900 and 
2300 GeV for 25 ftr 1 at 8 TeV, 300 fb” 1 and 3000 ftr 1 at 
14 TeV, respectively. We notice here that, for the general 
supersymmetric spectra realised in the pMSSM, values of 
the gluino mass below 1 TeV are still viable, and that the 
next runs will be able to extend the sensitivity beyond 
2 TeV. 

Since the gravitino LSP implies different constraints 




































9 



FIG. 4: Distribution of the gluino mass from the pMSSM 
scans with LSP gravitino and NLSP neutralino showing all 
the accepted points (black line) and those not excluded with 
25 fb -1 of 7+8 TeV data (grey line), 300 fb -1 (dashed line) 
and 3000 fb -1 of 14 TeV data (dotted line). The vertical 
lines indicate the masses at which more than 95% of our scan 
points with gluino mass below those values are, or will be, 
excluded by the LHC data. 


from dark matter, it is interesting to contrast the param¬ 
eter space allowed in this scenario to that of the MSSM 
with neutralino LSP, after the LHC searches. The main 
difference between the two models arises from the dark 
matter relic density constraint, which severely restricts 
the MSSM parameter space with neutralino LSP and 
modifies the occurrence of neutralinos of different na¬ 
ture and the relations of the masses of the other SUSY 
particles to that of the lightest neutralino. In the neu¬ 
tralino LSP scenario, bino-like neutralinos are relatively 
disfavoured by loose relic density constraint, because in 
general they lead to too large relic density. But imposing 
the tight relic density constraint highlights the bino-like 
scenario, while the wino- and higgsino-like scenarios be¬ 
come less attractive at the LHC since in these two cases 
the correct Vt x h 2 is naturally reached only for neutralino 
masses well beyond 1 TeV, thus pushing the SUSY spec¬ 
trum to high masses. In the wino scenario more com¬ 
plex decay chains of strongly-interacting SUSY particles 
occur also more frequently. In the bino scenario, it is 
necessary to have coannihilations, which means SUSY 
particles with masses close to that of the neutralino LSP 
are required. This has important consequences for the 
SUSY detectability at LHC, since compressed spectra or 
SUSY particles at high mass scales become favoured by 
the tight relic density bound in the neutralino DM sce¬ 
nario. 

Here, we use two sets of fl x h 2 constraints for the 
pMSSM points with neutralino LSP. First, we apply a 
tight bound requiring the two values to be in agree¬ 
ment allowing for systematic uncertainties, i.e. 0.090 < 
fl x h 2 < 0.163, and assuming that the LSP neutralino 


1 GravHIno LSP pMSSM~~| | Gravitino LSP cMSSM ~| 



Higgsino Higgsino 


FIG. 5: Nature of the lightest neutralino for pMSSM (top 
left) and cMSSM (top right) points with gravitino LSP and 
pMSSM with neutralino LSP with tight (bottom left) and 
loose (bottom right) relic DM constraints. 


saturates the observed dark matter. We also consider 
a looser bound by simply requesting that the neutralino 
relic density does not exceed the upper bound on the dark 
matter density from the PLANCK CMB, i.e. 10 -5 < 
fl x h 2 < 0.163, again after accounting for systematic un¬ 
certainties. This allows for other sources of dark matter 
in addition to the neutralino. Of the pMSSM points we 
have studied fulfilling the flavour physics and low-energy 
data constraints, 1.7% and 57% satisfy the tight and 
loose constraints, respectively. It is therefore clear that 
some amount of tuning is needed to obtain the right neu¬ 
tralino Dark Matter density. Moreover, the tight Q x h 2 
constraint gives only 6.5% of wino neutralinos, with 56% 
higgsinos and 26% bino. The loose constraint selects a 
sample of points with the lightest neutralino being wino 
and higgsino in 61% and 35% of the cases, respectively. 
This should be contrasted with the case of gravitino LSP 
models, characterised by an even distribution of points 
with different neutralino nature, from bino (38%) to wino 
(39%) and higgsino (22%) (see Figure [5]). In contrast 
to all these pMSSM scenarios, the gravitino model in 
the constrained MSSM is basically restricted to the case 
of bino-like neutralinos, which acccounts for 97% of the 
cMSSM accepted points, as shown in Figure [5] 


IV. RESULTS 

The viable regions of the MSSM parameter space with 
gravitino LSP excluded by the LHC Run-1 searches and 
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those projected for 300 fb -1 and 3000 fb -1 of data at 
14 TeV, in the case no excess of events above the SM 
backgrounds will be observed, have been studied for var¬ 
ious combinations of parameters. We present the results 
by computing the fractions of our scan points which are 


(will be) excluded by the LHC searches given in Table IV 


Quantitatively, the values of these fractions depend on 
the range of parameters used in the scans. However, qual¬ 
itatively the results hold when varying these ranges. The 
regions highlighted as being deeply probed by the LHC 
searches remain so even if the SUSY particle masses are 
pushed to higher values, although the fraction of points 
excluded will vary accordingly. A first observation de¬ 
rived from these results is the existence of a significant 
parameter space with neutralino NLSP when studying 
the pMSSM in sharp contrast with the results reported 
in previous studies using the cMSSM and fixed values of 
tan^ and tUq. The sensitivity of the LHC data in fal¬ 
sifying the currently viable MSSM points with gravitino 
LSP can be contrasted with that for models with neu- 
tralino LSP. In the previous section, we have analysed 
the difference in the nature of the lightest neutralino for 
the two models. Here, we concentrate on the fraction of 
points with SUSY masses up to 5 TeV excluded by the 
LHC searches. Figure [6] and Table [V| summarise our find¬ 
ings, with the regions in grey being not accessible with 
the statistics of our scans. The distributions of the vi¬ 
able points from the scans and their fraction excluded by 
a negative result with 300 fb -1 at 14 TeV in the grav¬ 
itino LSP and neutralino LSP scenarios with tight dark 
matter relic density constraint are shown by the figure in 
the plane of the gluino mass versus the neutralino mass. 
We see here that the reach in neutralino mass is simi¬ 
lar in the two cases and reaches approximately 1 TeV, 
while in the direction of the gluino mass, the LHC will 
be able to test better the large mass regions in the case of 
a neutralino LSP and DM. We trace this back to the fact 
that in case of a light neutralino DM, more supersymmet¬ 
ric particles have to be relatively light (e.g. the charged 
higgsino/winos for the case of higgsino/wino neutralino) 
and can give rise to detectable signals even if the coloured 
states are very heavy. 


The fractions of points excluded by the analyses at 8, 
14 TeV LHC and HL-LHC are given in Table |V| for grav¬ 
itino LSP and neutralino LSP with the tight and loose 
dark matter relic density constraints. Since the compo¬ 
sition of the sample of accepted points has significant 
differences in the nature of the lightest neutralino, as 
discussed above, these fractions are also given restricting 
the analysis to the points which have bino-, wino- and 
higgsino-like neutralino. We observe that the fraction of 
excluded gravitino LSP points is significantly larger (by 
^22% for the Run-1 data and almost 75% at 14 TeV) 
than that of models with neutralino LSP and tight relic 
dark matter constraints. This is largely due to the ef¬ 
fect of the higgsino-like neutralino component in the neu- 
tralino LSP sample. The higgsino/wino relic density nat¬ 
urally reaches the CMB value for neutralino LSP masses 



M(x) (GeV) 



M(x) (GeV) 

FIG. 6: Projected fractions of model points excluded by the 
LHC SUSY searches with 300 fb -1 at 14 TeV for the gravitino 
LSP (top) and the neutralino LSP scenario (bottom) in the 
plane defined by the gluino and lightest neutralino masses. 
For neutralino LSP scenario only points consistent with the 
tight U x h 2 constraints and with neutralino scattering cross 
section cr X p consistent with the LUX data are considered. The 
regions in grey have no points found in our scans. 


in the range 1.5-2.5 TeV. This pushes the SUSY spectrum 
to masses of ^2 TeV or more, i.e. at the limit of the LHC 
sensitivity or beyond. Therefore, if the Q x h 2 constraints 
are enforced, the LHC sensitivity to higgsino- and wino- 
like neutralinos LSP scenarios collapses. Indeed no points 
are excluded in this case from Run-1 data, as seen in Ta¬ 
ble 0 and also the future runs will not be able to cover 
the whole parameter space. Of course in case of the loose 
DM constraint, the neutralino mass is not pushed to such 
high values and the LHC is more effective in testing the 
model, at a level similar to that observed for the grav¬ 
itino LSP scenario. The LHC searches are also partic¬ 
ularly sensitive to the photino points, as well as to the 
other bino-like solutions. The fractions of these points 
removed by the Run-1 data and the projections to 300 
and 3000 fb -1 of data at 14 TeV are 37, 92 and 98%, 
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G LSP 

x° LSP 

.09 < Q x h 2 < .163 

x° LSP 

10“ 5 < n x h 2 < .163 

8 TeV 25 fb” 1 

All 

0.218 

0.100 

0.188 

Bino 

0.273 

0.268 

0.312 

Wino 

0.144 

0 

0.145 

Higgsino 

0.245 

0 

0.249 

14 TeV 300 ftr 1 
All 

0.745 

0.533 

0.694 

Bino 

0.835 

0.851 

0.864 

Wino 

0.614 

0.035 

0.615 

Higgsino 

0.808 

0.343 

0.811 

14 TeV 3 ab -1 
All 

0.845 

0.745 

0.806 

Bino 

0.917 

0.956 

0.927 

Wino 

0.733 

0.212 

0.736 

Higgsino 

0.902 

0.631 

0.907 


TABLE V: Fraction of pMSSM points with gravitino LSP and 
neutralino LSP fulfilling the flavour physics and low energy 
data constraints excluded by LHC searches. 


respectively. In the case of the gravitino scenario in the 
cMSSM, the projection for 300 fb -1 at 14 TeV has ^34% 
of the points excluded if no signal is observed, when re¬ 
stricting the analysis to points with SUSY particle masses 
below 5 TeV to cover a parameter space comparable to 
that used for the pMSSM. The difference in the fraction 
of excluded points is related to the distribution of the 
SUSY masses after applying all our selection criteria. 

The gravitino and neutralino LSP MSSM scenarios can 
be probed and differentiated by combining the LHC and 
DM direct detection searches. In fact, in the case of grav¬ 
itino LSP no signal from WIMP scattering on nucleons is 
expected in the underground DM searches. On the con¬ 
trary, if the neutralino LSP is the WIMP DM particle, 
a large fraction of the points to which ATLAS and CMS 
will become sensitive at the LHC forthcoming runs are 
also expected to give a signal at the LUX [6] or LZ DM 
experiments mm- The distribution of the fraction 
of points excluded by LHC with 300 fb -1 of 14 TeV data 
in the neutralino scattering cross section vs. neutralino 
mass plane is shown in Figure [7] for the gravitino and neu¬ 
tralino LSP and DM scenarios. We see here that, in the 
gravitino DM scenario with an NLSP mass below 1 TeV, 
a large part of the parameter points tested at LHC will 
also be tested by LZ and in that case the presence of a 
positive signal at collider without a direct detection coun¬ 
terpart will be a strong indication that the neutralino is 
not the LSP. For the case of a neutralino LSP and DM 
instead, the direct detection experiments are able to test 
also regions outside of the LHC sensitivity and detect 
heavy neutralinos, well beyond 1 TeV. Therefore, we ob¬ 
serve that only a small fraction of points in our scan are 
able to evade both the LHC and the direct detection con¬ 
straints. In quantitative terms, we find that 64 (62)% of 
the gravitino LSP and 82 (80)% of the neutralino LSP 
pMSSM points to which the LHC experiments will be 
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FIG. 7: Fractions of gravitino (top) and neutralino (bottom) 
LSP model points excluded by the LHC SUSY searches for 
the projection for 14 TeV with 300 fb -1 in the plane defined 
by the spin-independent neutralino scattering cross section 
and the neutralino mass. For the neutralino LSP the tight 
dark matter constraint is applied. 


sensitive after integrating 300 fb -1 (3000 fb -1 ) of data 
at 14 TeV could be also tested by the planned LZ direct 
detection experiment, when the loose Q x h 2 constraints 
are applied. In this case, it may be appropriate to rescale 
the x scattering cross section for the neutralino LSP case 
by the ratio of the neutralino relic density to the CMB 
value. This lowers the fraction of points given above to 
42 and 39 %, for 300 and 3000 fb -1 respectively, in line 
with the result that we obtain when enforcing the tight 
relic density constraints. 

Turning to the sensitivity over specific regions of the 
parameter space, the first set of variables we analyse is 
the relation of the neutralino relic density, mass and life¬ 
time (see Figures [8] and [9]). The sensitivity to the neu¬ 
tralino nature and mass is mapped, in part, onto these 
distributions. For example, the regions where the lightest 
neutralino is bino-like are preferentially excluded result¬ 
ing in the largest sensitivity from the LHC Run-1 data 
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FIG. 8: Fractions of gravitino LSP model points excluded by 
the LHC SUSY searches with the 7+8 TeV data (top) and the 
projection for 14 TeV with 300 fb -1 (centre) and 3000 fb -1 
(bottom) in the plane defined by the neutralino relic density 
and the neutralino lifetime. The BBN constraints obtained 
following ref. m are shown by the continuous lines corre¬ 
sponding to limits for neutralino hadronic branching fractions 
of 100%, 10%, 1% and 0% from bottom to top. 



M(xJ) (GeV) 


FIG. 9: Fractions of gravitino LSP model points excluded by 
the LHC SUSY searches with the 7+8 TeV data (top) and the 
projection for 14 TeV with 300 fb -1 (centre) and 3000 fb -1 
(bottom) in the plane defined by the neutralino mass and the 
neutralino lifetime. 


and the projections for the 14 TeV data being obtained 
for the parameter regions closest to the cosmological con¬ 
straints in Figure [8] Note on the other hand, comparing 
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with Figure [2j that also the region of wino NLSP will be 
well-tested in the future runs. Moreover, the points with 
larger lifetime, just at the boundary of the constraint, 
correspond on average also to smaller neutralino mass 
and are more easily tested at a collider. Such a trend is 
clearly visible in Fig. [9j 

In order to evaluate the sensitivity to the different as¬ 
sumptions used to derive the BBN constraints, we com¬ 
pare the results obtained using the constraints of ref. m 
which are used throughout our study, to those of ref. jl9j . 
We find that the differences are minimal and restricted 
only to very particular points at the boundary. So the 
fraction of pMSSM points with gravitino LSP fulfilling 
the flavour physics and low energy data constraints ex¬ 
cluded by LHC searches for these two implementations 
of the BBN bounds agree to better than 5%. The two 
samples of points differ mostly in the particular value of 
the maximal gravitino mass, which does not affect the 
LHC phenomenology. Moreover, since we are not scan¬ 
ning systematically on that mass, but taking a compara¬ 
ble number of points at the boundary, such a difference 
does not influence substantially the fraction of excluded 
points. 

The fractions of excluded points in the plane of the 
hadronic branching fraction as a function of the neu- 
tralino mass are shown in Figure [ 10 | Here, it is interest¬ 
ing to observe that the photino solutions can be largely 
excluded in the case the HL-LHC will not observe an ex¬ 
cess of events in the SUSY searches: 97% of the points 
in our scan are excluded with 3000 fb -1 at 14 TeV, as 
shown in the right panel of Figure [lO] But in general the 
sensitivity of the LHC in the gravitino DM scenario is 
again mostly determined by the neutralino NLSP mass 
scale rather than by its composition. 

The LHC is also particularly sensitive to the test of 
gravitino LSP scenarios with relatively light NLSP neu- 
tralinos and large reheating temperature, T##, close to 
the constraints derived from thermal leptogenesis, for 
which we assume the lower bounds of 2 x 10 9 GeV [381 HD] 
or 5 x 10 8 GeV, in presence of flavour effects [iHIST] , as 
shown in Figure [Tl] Compared to cMSSM studies, which 
have essentially only bino x NLSP as we have discussed 
above, the pMSSM scans generated for this analysis have 
points with wino x LSP that enable to increase the values 
of Trh reachable in this scenario (see Figure |3| and make 
the LHC sensitivity compelling for their test. These re¬ 
sults are consistent with those obtained by the analysis 
of Ref. [73] that found a reheating temperature of order 
of 10 9 GeV for the case of 1.25 TeV gluino and maxi¬ 
mal gravitino mass around 70 GeV. In our scan we also 
reach larger gravitino masses, i.e. also larger Trh , since 
the neutralino relic density in the pMSSM reaches lower 
values than those obtained in m- 

Indeed a high reheating temperature is necessary to ob¬ 
tain the right gravitino abundance when the neutralino 
is light and has too low relic density to produce a sizeable 
gravitino population in its decay. We see clearly in the 
figure also the complementary region of parameter space, 





M(xJ) (GeV) 


FIG. 10: Fractions of gravitino LSP points excluded by the 
LHC SUSY searches with the 7+8 TeV data (top) and the 
projection for 14 TeV with 300 fb -1 (centre) and 3000 fb -1 
(bottom) in the plane defined by the neutralino hadronic de¬ 
cay branching fraction and the neutralino mass. Note the 
photino solutions at small branching fraction values being al¬ 
most completely excluded in case of no SUSY signal at the 
HL-LHC. 
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FIG. 11: Fractions of gravitino LSP points excluded by the 
LHC SUSY searches with the 7+8 TeV data (top) and the 
projection for 14 TeV with 300 fb -1 (centre) and 3000 fb -1 
(bottom) in the plane defined by the reheating temperature 
and the neutralino mass. The horizontal lines show the re¬ 
heating temperature constraints from leptogenesis. 


where the whole gravitino density is indeed generated by 
the light neutralino decay and therefore the reheating 


temperature has to be very low. Both these two param¬ 
eter regions will be completely tested in the next run of 
the LHC for neutralino masses below 500-700 GeV. 

In these scenarios, the LHC sensitivity mostly comes 
from the constraints on the gluinos which have impor¬ 
tant implications on the viable values of the reheating 
temperature. The interplay of the gravitino and gluino 
masses and Trh for the points of our scan after imposing 
the current and expected LHC constraints in absence of 
a signal is highlighted in Figures [I2| and 

From Fig. [l2j we can see clearly that the requirement of 
a specific reheating temperature and the LHC exclusion 
on the gluino mass are completely orthogonal constraints. 
In the upper Trh region we point out that the bound¬ 
ary of the sample points is determined by the overclosure 
bound QqH 2 < 0.163: since the gravitino density is in¬ 
creased both by large reheating temperature and gluino 
mass, as given in eq. ([l]) , such constraint imposes an an¬ 
ticorrelation between the maximal temperature allowed 
and the gluino mass. Note that the fact that the low 
reheating temperature region is tested also at very large 
gluino masses, well beyond the LHC reach, is connected 
to the presence of a parameter region where the graviti- 
nos are produced via neutralino decay and in that case 
neither the LHC constraints nor the gravitino density are 
derived from the coloured sector of the model. In Fig. [13] 
we can see instead the points in our scan along the grav¬ 
itino mass. Since we took mostly values of the gravitino 
mass around the BBN bounds, the distribution is not 
uniform in this direction. We note that here the points 
corresponding to low Trh have mostly a very small grav¬ 
itino mass. The particular behaviour of the upper bound 
is not due, in this case, directly to the overclosure bound, 
which would give T^ x oc M^, but to the interplay of the 
overclosure bound with the BBN constraints. 

As we have discussed, the contribution of heavy gluinos 
to the production of gravitinos after inflation affects the 
QqH 2 value, as can be seen from eq. |lj) and LHC bounds, 
if combined with the additional requirement of success¬ 
ful thermal leptogenesis may effectively test the result¬ 
ing gravitino LSP scenarios [32, 39;. This is highlighted 
in Figure [l4j which summarises the results of our scans 
by giving the contours for the minimum values of the 
gravitino masses obtained as a function of the gluino 
mass for different choices of Trh- This plot can be 
read in two complementary ways. First, by requiring a 
particular minimal reheating temperature together with 
the relic density constraint, we can set strong bounds 
on the gluino and gravitino masses, since these deter¬ 
mine the gravitino abundance in most of our parameter 
space. For gluinos heavier than 800 GeV, satisfying the 
present LHC exclusion bounds, the gravitino mass has 
also to be substantial for reheating temperature above 
10 9 GeV. Conversely, for each gravitino mass value, a 
choice of Trh corresponds to an upper bound on the 
gluino mass to avoid overclosure. Then, by comparing 
these contours with the reheating temperature needed 
by thermal leptogenesis discussed above, we observe that 
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FIG. 12: Fractions of gravitino LSP points excluded by the 
LHC SUSY searches with the 7+8 TeV data (top) and the 
projection for 14 TeV with 300 fb -1 (centre) and 3000 fb -1 
(bottom) in the plane defined by the reheating temperature 
and the gluino mass. The horizontal lines show the reheat¬ 
ing temperature constraints from leptogenesis and the vertical 
lines give the gluino mass value for which more than 95% of 
the points below it are excluded by the LHC SUSY searches. 


FIG. 13: Fractions of gravitino LSP points excluded by the 
LHC SUSY searches with the 7+8 TeV data (top) and the 
projection for 14 TeV with 300 fb -1 (centre) and 3000 fb -1 
(bottom) in the plane defined by the reheating temperature 
and the gravitino mass. The horizontal lines show the reheat¬ 
ing temperature constraints from leptogenesis. 


the LHC data should either observe a gluino signal or 
exclude the gravitino DM scenarios with gravitino mass 
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FIG. 14: Minimum values of the gravitino mass obtained in 
the gravitino LSP scans as a function of the gluino mass for 
different choices of the reheating temperature (in GeV). The 
values of 2 x 10 9 and 5 X 10 8 GeV correspond to the limits 
imposed by leptogenesis |3HI EDI HUES- The thin diagonal 
line indicates the lower limit of the region where the gravitino 
is not the LSP. The vertical lines give the values of the gluino 
masses for which more than 95% of the scan points below 
them are excluded by the LHC SUSY searches at 7+8 TeV 
(continuous) and projected at 14 TeV for 300 (dashed) and 
3000 fb -1 (dotted), if no signal is observed. 


below 200-300 GeV and classical thermal leptogenesis, 
after accumulating sufficient statistics. In this study the 
availability of the data set anticipated from the HL-LHC 
operation is of crucial importance, since the sensitivity 
afforded by 3 ab -1 of data at 14 TeV would cover the 
region compatible with the leptogenesis data up to at 
least Mq ~ 300 GeV for a reheating temperature bound 
of Trh > 2 x 10 9 GeV and more for higher values of 
Trr. In this case, the lack of observation of a gluino in 
the LHC data would basically exclude most of the grav¬ 
itino LSP scenarios discussed here for the reheating tem¬ 
peratures required by thermal leptogenesis. Therefore, 
gravitino LSP models, in conjunction with constraints 
derived from leptogenesis, offer a strong motivation for 
the HL-LHC program. Scenarios of flavoured leptoge¬ 
nesis, allowing for lower Trh > 5 x 10 8 GeV, may still 
accommodate gluino masses even larger than those which 
can be probed at the HL-LHC, but could still be probed 
for gravitino masses below ~ 60 GeV. 


V. CONCLUSIONS 

The search for dark matter is tightly connected to that 
of new physics beyond the SM. These searches bring 
together the efforts and the data accumulated by the 
LHC experiments but also the DM underground direct 
searches. Gravitino Dark Matter represents a compelling 
scenario in Supersymmetry. In this paper we have pre¬ 
sented a detailed study of the gravitino LSP and Dark 


Matter scenario within the phenomenological MSSM and 
contrasted it with the solutions available in the CMSSM. 
Cosmological bounds on the scenario result in upper 
bounds on the gravitino mass and the reheating temper¬ 
ature necessary to fulfil the requirement of a DM energy 
density in agreement with CMB data. 


Our analysis focused on the large gravitino mass re¬ 
gion, at the boundary of the nucleosynthesis constraints, 
because that part of the parameter space allows for the 
largest possible values of the reheating temperature. The 
present and projected LHC constraints considered in this 
study include the monojet bounds that are important to 
test the parameter regions with large mass degeneracy, 
in particular the case of higgsino and wino NLSP. The 
use of the unconstrained pMSSM model with 19+1 free 
parameters has scenarios with bino-, wino- and higgsino- 
like NLSP, which are not available to studies carried out 
in the constrained MSSM with 5 parameters which yield 
almost uniquely solutions with bino-like NLSP neutrali- 
nos. We find that the scenario with gravitino LSP dark 
matter in the pMSSM is characterised by a very differ¬ 
ent composition of the neutralino NLSP in comparison to 
the case of a neutralino LSP. In particular all three types 
of compositions are nearly equivalently possible in the 
gravitino DM scenario, while for a WIMP neutralino a 
higgsino composition is preferred. Such a specific compo¬ 
sition is unfortunately more difficult to test at a collider, 
since the DM mass scale tends to be above 1 TeV. On the 
other hand, direct detection experiments could give a fu¬ 
ture detection, even for large neutralino mass, and cover 
a large part of the yet allowed parameter space. But if 
no signal appears in searches at those experiments and 
in direct production at the LHC, strong constraints on 
the neutralino WIMP scenario will be obtained, making 
the gravitino LSP scenarios particularly interesting. We 
have seen that in case of gravitino DM and neutralino 
NLSP, the absence of a direct detection signal should al¬ 
low to exclude the neutralino DM scenario and point to 
an unstable NLSP in a large part of the parameter space. 


The interplay of the LHC and cosmology data can be 
pursued even further in gravitino dark matter scenarios. 
The relation between the gravitino relic density and the 
gluino or gaugino masses, provides us with an efficient 
way to test also the cosmological scenario. The region of 
the pMSSM parameter space with large gluino masses, 
beyond the projected LHC sensitivity, corresponds to 
gravitino energy densities that are too large at the high 
reheating temperatures required by thermal leptogenesis. 
Therefore, if the next runs of the LHC and especially the 
High-Luminosity operation will not observe any signal 
of gluino direct production, the high reheating tempera¬ 
ture region will be excluded for gravitino masses below 
200 GeV. 
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